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Abstract 


A series of solid solution «TiFe,O,.+(1—.«)Fe,0, was synthesized by ceramic 
method over a whole range of the composition, 1=x=0. Changes of the Curie 
point, saturation moment and the lattice parameter with the composition were 
examined. Generalized titanomagnetite having some vacant site in the structure 
normally occupied in a spinel was also prepared by oxidizing the TiFe,O,-Fe,0, 
solid solution series. The region within which the spinel structure can be in 
existence as a single phase was settled in a FeO—Fe,0,-TiO, system. An equal 
lattice parameter line, equal Curie point line and equal saturation moment line 


for the generalized titanomagnetite were drawn on the FeOQ-Fe,0,-TiO, system. 


1. Introduction 


It has been established from recent study on the ferromagnetic minerals contained 
in rocks that the chemical composition of the natural titanomagnetite with spinel 
structure does not always accord with TiFe,O,-Fe,O, join but in most cases deviates 
from this join towards TiFeO,-Fe,O; join side and occasionally FeO side in a FeO- 
Fe,0,-TiO, system. As long as we deal with these natural specimens, it may be very 
probable that the various magnetic properties are considerably scattered with respect 
to the calculated TiFe,O, content which was tentatively used in the previous paper as 
a measure representative of their chemical composition (Akimoto; 1954, 1955). 

It is indispensable for interpreting correctly the magnetic properties of these 
natural titanomagnetites to perform a systematic study on the synthetic titanomagnetites 
which were prepared under the known physico-chemical condition. Several magneto- 
chemical studies along this line have already been reported by Ernst (1943), Pouillard 
- (1950), Chevallier and Girard (1950), Kawai, Kume and Sasajima (1954), and Chevallier, 
Bolfa and Mathieu (1955). The data obtained hitherto by these authors are still 
unsatisfactory for applying to the interpretation of the magnetic properties of natural 
titanomagnetites. Pouillard’s synthetic titanomagnetites are limited to the solid solution 
xTiFe,O,.-(1—x«)Fe,0, of TiFe,O, and Fe,;O, over a range of 0.42>x>0. He synthesized 
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his specimens either by heating mixture of TiO, and Fe,O, in an evacuated furnace 
for one hour at 900°C or by reducing mixture of TiO, and Fe,O; with hydrogen gas 
for two hours at 450°C. In each case other phases were also formed and he estimated 
the molecular percentages of TiFe,O, by applying the law of Vegard (linear relation 
between crystal parameter and molecular camposition). Therefore, as Nicholls (1955) 
pointed out in his review on rock-forming ferromagnetic minerals, there remain still 
some ambiguities upon the chemical composition on his specimens. Although Kawai, 
Kume and Sasajima have reported that they succeeded in extending the range of synthetic 
spinels to x=0.8 with the aid of usual ceramic method, they do not give any detailed 
description on the way in which the composition of the synthetic specimen was 
determined. 

On the other hand, very little work has been done with respect to the magneto- 
chemistry on the cubic solid solutions between TiFe,O,-Fe,O, join and TiFeO;-Fe,O; 
join. The only information available is the work carried out on TiFeO,;-Fe;O, solid 
solution series by Chevallier and Girard (1950) and Chevallier, Bolfa and Mathieu (1955). 
They have synthesized cubic solid solutions intermediate between Fe,O, and TiFeO; 
up to 37 molecular percent of TiFeO; by means of a borax method. 

Under the circumstances, we intended to extend the studies hitherto made to the 
more broad region between TiFe,O,-Fe,0, join and TiFeO,-Fe,O; join. As the first 
step of the study we synthesized the solid solution xTiFe,O,-(1—x)Fe;0, over a whole 
range of the composition. Then we synthesized the cubic solid solution in the concerned 
region by oxidizing the TiFe,O,-Fe,O, series specimens. The particular interests shown 
in this paper are related to the following two: 

(1) The one is to settle the region on the FeO-Fe,0;-TiO, system, within which 
the spinel structure can be in existence as a single phase. 

(2) The other is to draw an equal lattice parameter line, equal Curie point line 
or equal saturation moment line on the FeO-Fe,O;-TiO, ternary system. 


2. Preparation of TiFe,0,-Fe,0, Solid Solution Series 


The specimens dealt with in the present study were prepared by the sintering 
procedure generally used in the ceramics. A fine powder mixture of pure Fe,O,, TiO, and 
electrolytic iron in the desired proportion, sealed in a silica tube evacuated approximately 
to 10°*mmHg, was quenched from 1150°C after being maintained at this temperature 
for six hours. 

Twenty specimens were prepared by this method over the whole range of 
1=>x=0. The chemical composition of the synthesized specimen is plotted by full circles 
on the FeO-Fe,0;-TiO, diagram in Fig. 1 and is given in Table I. 

Analysis by means of the “Norelco” X-ray diffractometer shows that these 
specimens are generally. composed of a single phase having the spinel structure. No 
trace of the ilmenite-hematite series could be found except the specimen S 2-4, HICH: 
in this specimen the amounts of the ilmenite-hematite series minerals were estimated 
to be less than ten percent from the careful examination of the diffraction chart, -_; 
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Lattice parameters of all synthetic specimens determined from the X-ray analysis are 
also given in Table I. In Fig. 2 the lattice parameter is plotted against the TiFe,O, 
content x, which was tentatively determined from the data of chemical analysis as a 
mean value of the following two values of x derived from different assumptions: the 
one is based on the assumption that all the titanium oxides are involved into the 
ulvospinel (TiFe,O,), the other being based on the assumption that all the ferric oxides 
are attributed to magnetite (Fe,O,). As will be seen in Fig. 2 and Table I, it was 
established that a solid solution xTiFe,O,+(1—x)Fe,;O, of spinel structure can be formed 
continuously throughout the whole range of 1=x=0. The lattice parameter varies 
almost linearly with the composition from about 8.39 A of Fe,0, to 8.53 A of TiFe,O,. 


3. Oxidation of TiFe,O,-Fe,0, Solid Solution Series 


Some experimental work on the oxidation of titanomagnetite has been carried 
out by Nagata and Ozima (1955) in connection with the interpretation of a particular 
phenomenon of thermoremanent magnetism, i.e. anomalous increase of thermoremanent 
magnetism. They suggested that a material having spinel structure could be remained 
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Fig. 2 Relation between lattice parameter and chemical composition 
in TiFe,O,-Fe;O, solid solution series. 
as a single phase in the region between TiFe,O,-Fe;O, join and TiFeO,-Fe,O, join. 
We also obtained the cubic solid solutions in the concerned region by oxidizing the 
TiFe,O,-Fe,O0, series specimen. 

The specimens of TiFe,O,-Fe;O0, solid solution series, packed in an open air 
porcelain boat, were inserted into an electric furnace after the temperature in the 
furnace had settled at the designated temperature, and after being maintained there 
for a designated length of time, taken out and cooled in air. The length of the heat 
treatment was two hours in most cases and the temperature of heat treatment was the 

intermediate between 400°C and 550°C. 

Ten specimens from S 2-1 to S2-10 and other four specimens S 5-1, $ 5-3, S5-5 
and S5-6, described in just preceding section, were used for the present purpose. 
Chemical compositions of the heat-treated specimens are also plotted by hollow circles 
in Fig. 1 and given in Table II. The heat-treated specimens are distinguished in the 
figure or table by a mark ’or ” from the original ones. As seen in Fig. 1, the chemical 
composition of the original specimen was shifted by the heat treatment towards the 
directions of arrows presented in the figure along the theoretical reduction-oxidation 
line. 

The X-ray analysis by the “Norelco” X-ray diffractometer revealed that the 
specimens still kept their spinel structure after the heat treatment. As for the greater 
parts of the specimens (S2-3’, $2-5’, S2-6’, S2-7/, S2-8’, S2-9’, S5-1’, S5-3’, S5-5’, S5-6/, 

S5-1”, S5-3”, S5-5” and S5-6”) no trace of the ilmenite-hematite series minerals could 
be detected on the diffraction chart. But as for the specimen S2-1’, fairly large 
amounts of ilmenite-hematite series minerals were found. The amounts were estimated 
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from these experimental results that the spinel phase can be in existence as a single 
phase in fairly broad region between the TiFe,O,-Fe;O, join and the TiFeO,-Fe,O; join 
under an appropriate condition of temperature and oxygen pressure. It may be very 
probable to expect that the region of the spinel phase can be extended more close to 
the TiFeO,;-Fe,O; join if we find the appropriate heat treatment. The fact that a 
complete solid solution can be prepared between Fe,O, and 7 Fe,O;, (Hagg; 1935) may 
also support the above possibility. 


4. Thermomagnetic Curve and Curie Temperature 


The variation in magnetic moment of all the synthesized spinel specimens as a 
function of temperature in a constant magnetic field of a few thousands Oersteds was 
measured by means of a magnetic balance described in the previous paper (Akimoto; 
1954). In the present study the measurement was carried out under the evacuated 
state of about 10°*mmHg. The Curie temperature determined from the thermomagnetic 
curve is listed in Table I and Table II. As for the specimens of the TiFe,O,-Fe,0, 
solid solution series, the Curie temperature is plotted against chemical composition in 
Fig. 3, where it can be found that the Curie temperature changes almost linearly with 
composition from 580°C of Fe,0, down to about —150°C of TiFe,Q,. 

As for the specimens with lower Curie temperature the measurement of the 
thermomagnetic curve down to the liquid nitrogen temperature in a field strength of 
8450 Oe was practised by the help of Ishikawa at Institute of Science and Technology, 
Tokyo University. A few examples of the thermomagnetic curve are shown in Fig. 4, 
where we can see the peculiar mode of thermomagnetic curve like Néel’s P-type for 
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Fig. 4 A few examples of the thermomagnetic curve of TiFe,O,-Fe,O, 

solid solution series. 


the specimens S2-8 and S2-9. This may be a remarkable contrast to the fact that the 
thermomagnetic curve of the specimens of which TiFe,O, content is less than 60 percent 
does not show any peculiar mode but is very similar to that of Fe,O,. It must be 
also noticed from Fig. 4 that the specimen of which chemical composition is determined 
to be 0.96 TiFe,O,-0.04 Fe,;0,(S1-4) becomes ferrimagnetic below 170°K, its intensity of 
magnetization at 80°K in 8450 Oe being about 4 e.m.u./gr. This suggests that pure 
TiFe,O, becomes ferrimagnetic below about 120°K. So long as the inverse spinel 
structure is kept in TiFe,O,: that is, titanium ions are confined to the octahedral sites 


in the spinel structure, the intensity of saturation magnetization ought to be zero at =~ 


0°K. The present situation that TiFe,O, becomes ferrimagnetic at the lower tempera- 
tures may suggest that some titanium ions enter into the tetrahedral sites with the 
result that the compensation of magnetic moment between tetrahedral and octahedral 
sites is disturbed. Such disordering of cations from their theoretical distribution is 
very pr apabie “Ms oe eee TiFe,0,-Fe,0, series” 
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field strength of about 3000 Oe. The intensity of saturation magnetization was deter- 
mined by extrapolating the following empirical formula to H= oo 
o=o0(t—2), 
H 

where A is the constant determined experimentally for every specimen. As for the 
specimens with the lower Curie temperature the magnetization curve was measured at 
the liquid nitrogen temperature by the help of Ishikawa. The intensity of saturation 
magnetization of the synthesized titanomagnetite at room temperature or 80°K is given 
in Table I and Table II. The intensity of saturation magnetization of these specimens 
at 0°K was estimated by extrapolating their thermomagnetic curve to 0°K. These 
estimated saturation moments of the TiFe,O,-Fe;0, series specimens are also shown in 
Table I with the Bohr magneton number per one molecule of «TiFe,O,-(1—x)Fe,O,. 
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calculated (Akimoto; 1954, 1955). The other configuration of cation was postulated by 
Néel (1955) and Chevallier et al. (1955) by taking the Verwey’s empirical law concerning 


the location of cations in the spinel structure into consideration, that is 
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for which the saturation moments of 4—6x (for i<, <5] and 2—2x (for x>>) are 
calculated. 

The actual values of the saturation moments of the synthetic TiFe,O,-Fe,O, 
solid solution series are the intermediate between these two. It may also be worthwhile 
to note that the actual values in our synthetic specimens were always larger than the 
Néel-Chevallier’s theoretical value over the whole range of the composition. This seems 
to suggest that some titanium ions enter into the tetrahedral sites apart from the 


normal octahedral sites of Fe®+ (Fe2?t,Fe3t,,Tit*)O3- or Fe3+,,Fe3?_, (Fe3t,Tiz*)O7-. In 
these cases, as suggested by Gorter (1957), the cation distribution should be considered 
as 

Feit, Tig'(Fej,.Fejto, ce 1 i354), 


or Fejt, Vetinio, lig Mereegeieea loa, 


the saturation moment being 4—6x+10a@ and 2—2x+8a respectively. Then, the 
saturation moment should be increased with the increase of Ti ions in the tetrahedral 
sites. As for the present synthetic specimens prepared by the quenching method from 
high temperature, such disordering of cations may be very probable. Similar effect of 
the disordering of cations from the ideal inverse spinel structure has already been 
demonstrated for MgFe,O, by Bertaut (1952) and for NiFeAlO, by Gorter (1954). The 
implications of such a behaviour for rock magnetism are discussed by Néel (1955) and 
Verhoogen (1956) in connection with the mechanism causing self-reversal of the remanent 
magnetization of ferromagnetic minerals. The effect of annealing on the saturation 
moment of the synthetic specimens is now under investigation. 

The variation of saturation moment of the TiFe,0,—Fe,O, solid solution series 
with oxidation was also investigated. Comparing the saturation moment of the original 
specimen (Table I) with the oxidized one (Table II), we may notice the following 
remarkable situations. That is, the saturation moment at the room temperature of 
the specimen with the lower content of TiFe,O, decreases considerably according as 
the degree of oxidation proceeds, but that of the specimen with higher content of 
TiFe,O, is increased by oxidation. These experimental results would be explained fairly 
satisfactorily if we regard the specimen in the intermediate between TiFe,O,-Fe;,0, 
join and TiF eO;-Fe,O; join as a generalized titanomagnetite having some vacant sites 
in the structure normally occupied by metallic ions in a spinel and if we assume that 
the solid solution can be in existence between any composition on the. TiFe,O,-Fe,O, 
join and the corresponding oxidation product on the TiFeO,-Fe,O,; join. It is needless 
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to say that we must assume here the cubic crystal form on the TiFeO;-Fe,O; join. As 
for Fe,O3, y Fe,O; is well known and its magnetic structure is considered to be 


3+ 8+ 2- x = a a : : . : 
Fe (Fes 010% , Where [] denotes the vacant site in the structure normally occupied in a 


spinel. On the other hand as for the cubic form of TiFeO; magnetic structures of 


Fee" (Fes Ti *)OF and Fe?+(TittFe?*D 


Or have already been postulated by 


Chevallier et al. (1955) and Nicholls (1955) respectively for interpreting the saturation 
moment of the TiFeO;-Fe,O, solid solution series. In the present paper hypothetical 
cubic form of the TiFeO,-Fe,O; system, yTiFéO;-(1—y)Fe,O; is taken to be 


)0%- 


in making use of the Chevallier’s ;TiFeO;. This implies that the saturation moment 


Fe3+,Fe3* 8+ 2+ Ty 4+ 

: Fer, ry Pegi _y Fes, Tis 41-y) 
of the cubic TiFeO;-Fe,0; system changes linearly with the composition according to 
equation 


MMs 


We showed the variation of saturation moment along the reduction-oxidation line in 
Table III, where the original value of saturation moment of TiFe,O,-Fe,O, join and the 
final value of TiFeO;-Fe,O; join are given for different values of Fe/Fe+Ti in the 
unit of Bohr magneton number per molecule. As for the saturation moment of 


Table III Change in saturation moment of TiFe,O,-Fe,O, solid 
solution series with oxidation. 


; | Original spinel, Oxidation product, 
Fe/Fe+ Ti xTiFe,O,-(1—x)Fe,O,° yTiFeO;-(1— y)Fe,O3 
ee saturation moment (/tz) y saturation moment(/Zz) 
1.000 0 4.0* 0 S203 
0.966 0.1 3.4 (3.6)** 0.067 Sn EL 
0.933 0.2 Aste! (3:2) Ontss: 2.89 
0.900 0.3 Das (2.8) 0.200 2260 
0.866 0.4 1.6 (2.4) 0.267 2.44 
0.833 0.5 1.0 (2.0) Of333 Dery 
0.800 0.6 0.8 (1.6) 0.400 2.00 
0.766 0.7 0.6 (1.2) 0.467 1.78 
0.733 0.8 0.4 (0.8) 0.533 1.56 
0.700 0.9 0.2 (0.4) 0.600 ISS 
0.666 1.0 0 0.667 i abil 


* after Néel-Chevallier’s model. 
** after Akimoto’s model. 


TiFe,O,-Fe,O, join, two sets of values are listed in the table; ie. Néel-Chevallier’s 
value (Néel; 1955, Chevallier et al.; 1955) and the value postulated previously by one of 
the authors (Akimoto; 1954). From the table it is clearly seen that the saturation 
“moment of the specimen with lower content of TiFe,O, decreases according as the © 


oxidation proceeds, while that of the higher content of TiFe,O, increases. This situation 
accords well at least qualitatively with the above-mentioned experimental results, 
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6. Conclusive Remarks 


The accumulation of the reliable data on the chemical, crystallographic and 
magnetic properties of the spinel specimens enables us to draw an equal lattice constant 
line, equal Curie point line or equal saturation moment line on the spinel region in the 
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Fig. 8 Lines of equal saturation moment at the room temperature 

of the generalized titanomagnetite in a FeO-Fe,O,-TiO, system. 

FeO-Fe,0;-TiO; system. Figs. 6, 7 and 8 are the equal lattice parameter diagram, 
equal Curie point diagram and equal saturation moment diagram respectively. The 
results of Chevallier and Girard’s study (1950) are also taken into account in drawing 
these diagrams. ; 

It can be easily found from these diagrams that the general view of the equal 
lattice parameter diagram and the equal Curie point diagram resembles very much 
with each other; that is, the parallelism of the equal lattice parameter line and the 
equal Curie point line is generally found. Another remarkable character found in these 
diagrams is that the equal Curie point line representing the temperature higher than 
about 500°C has a marked curvature in contrast to the straight line of the lower 
temperatures. 

It must also be mentioned here that such an attempt to draw an equal Curie 
point line on the FeO-Fe,O0,-TiO, system has already been done by Chevallier and 
Girard (1950). They suggested the Curie point variation may be expressed by 

6°K (Curie point)=1245s 
0 
“eek ze. aor 
The present diagram shown in Fig. 7 accords well with that postulated by them within 
the Fe,0,-TiFeO,-TiFe,O, system except higher temperature side. 

The general view of the equal saturation moment diagram differs greatly from 
other two diagrams, that is, the equal saturation moment line intersects the equal 
lattice parameter line or the equal Curie point line in the FeO-Fe,O0;-TiO, diagram. 

These diagrams seem to have a deep significance especially in the field of 
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geological application of rock magnetism. It has already been shown that the magnetic 
analysis, together with chemical and X-ray analyses, is sometimes useful for deter- 
mining the constitution of the ferromagnetic minerals contained in rocks. As the 
present diagrams show, however, the chemical composition having a definite lattice 
parameter or Curie point is not represented by a point but form a line on the 
FeO-Fe,0;-TiO, system. Hence, even if the lattice parameter or Curie point of the 
ferromagnetic minerals were determined with sufficient accuracy by X-ray or magnetic 
analysis, it is impossible to fix their chemical composition to a point on the 
FeO-Fe,O;-TiO, system. The determination of the chemical composition by these method 
is always accompanied with some uncertainties. When we measure the saturation moment 
of the specimen at the same time, this uncertainty will vanish, the chemical composi- 
tion of the specimen being able to be fixed to a single point. It is also expected that 
the reliable data on the specific gravity of the specimen is useful for the determination 


of the chemical compositions. 


In concluding the present paper, the authors wish to express their sincere thanks 
to Prof. T. Nagata for his constant guidance throughout the study. They also wish 
to thank Prof. I. Iwasaki for his kind encouragement and interest. The authors are 
also indebted to Dr. Y. Ishikawa for his magnetic measurement at lower temperatures. 
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The Variation Type of the Atmospheric Potential Gradient 
and the Geographic latitude of the Globe 


By- Miyoji GOTO 


Aso Institute, Kyoto University 
(Read May 16, 1953; Received March 20, 1958) 


The global distribution of the annual and diurnal variation of the atmospheric 
potential gradient is discussed, using the material which is derived from H. Benndorf 
(1928) as shown in Table 1. Making use of A (Jahresamplitude des P.G. in V/m) and 


Table 1 
Lauf- Jahres- aay Ge 
Geogr. Beobach- | Mittel | tude | tude 
ende Beobachtungsort des. des des A/D 
Nr. Breite tungszeit vin ae Ee 
V/m(A)|V/m(D) 

it Ebeltofthafen 79.1°N 1913-1914 95 76 16 4.75 

& Karasjok 69.3°N | 1903-1904] 139 119 110 1.08 

3 Vassijaure 68.4°N | 1909-1910 89 110 50 BD 

4 Upsala 59.9°N | 1912-1914 70 59 50 1.18 

5 Aas 59.7°N | 1916-1923] 104 105° 46 2.28 

6 Eskdalemuir 55.3°N | 1914-1920} 263 161 109 ol 

fi Potsdam 52.4°N 1904-1923} 202 IWeAg Thee 1.76 

8 Kew 51.5°N | 1898-1920; 317 239 130 1.83 

9 Val Joyeux 48.8°N | 1923-1924 90 63 48 SSL 
10 Miinchen 48.1°N | 1905-1910] 168 161 132 1522 
11 Kremsmunster 48.1°N | 1902-1916) 105 79 60 Wee 
WZ; Davos 46.8°N | 1908-1910 64 68 44 1.54 
13 Triest 45.6°N_ | 1902-1905 73 25 86 0.29 
14 Tortosa 40.8°N | 1910-1924} 106 oe 57 0.65 
15 Washington 39.0°N | 1920-1923] 179 - — = 
16 Helwan 29.8°N | 1909-1914) 150 54 61 0.88 
17 Batavia 632° 1890-1900} 120 53 145 0.36 
18 Apia 1328 1913-1924) 112 22 92 0.24 
19 Rio de Janeiro 22.9°S | 1910-1914} 128 80 — = 
20 Buenos Aires 34.5°S | 1911-1912} 136 78 101 0.77 
21 Melbourne LY tc eae 1858-1860 145 15 125 0.60 
he Petermannsinsel 6522-5 1909 164 190 92 2.06 
23 Cap Evans 7156" o=- TOl1=1 912 ey 40 30 Ty33 
24 | Mc Murdo Sund 77.8°S | 1902-1903] 93 63 = 5 
25 | Auf dem Ozean 80°N | 1915-1921) 124 20 43 | 0.46 

(Carnegie fahrt) 60°S 


(179 ) 
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D (Tagesamplitude des P.G. in V/m) from the Table, the ratio A/D for each stations 
can be calculated numerically as indicated beside the Table 1. Then the whole stations 
on the globe can be divided into two groups according to A/D>1 and A/D<1. 

The relation between the ratio A/D and geographic latitude of each stations is 
shown graphically in Fig. 1. From Fig. 1*we can see distinctly that the stations in 
which A/D>1 are all located in the high latitude regions, northwards from 46°N and 
southwards from 46°S in both hemisphers, and the remaining stations in which A/D<1 
are all located in the low latitude zone from 46°N to 46°S. Judging from A/D=0.46<1 
from the data of the Carnegie cruise on the ocean from 80°N to 60°S, it is thought that 
these data may be mainly from the wide sea area of the low latitude zone.. As above 
mentioned, the high latitude caps and low latitude belt can be distinguished distinctly 
with respect to the annual and diurnal variations of the atmospheric potential gradient. 
It is desirable to make an observation in the neighbourhood of the transitional latitude 
(46°) in both hemispheres, especially in southern hemisphere. 


A/D 
2.5 


2.0 


tN og 


is 
*s ; o} ar | ; 


« ~ — a | 
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(Aquatorialgiirtel) indicate somewhat different mode of annual variation, maximum in 
Juli-Aug and minimum in Nov-Marz. 

It is supposed that the mode of the annual variation of the potential gradient 
suggested by H. Benndorf and the global distribution of the ratio A/D have deep 
mutual connection each other. it is too hasty to state the unique interpretation 
governing the above mentioned world wide phenomena, but we may say in moderation 
that judging from A/D>1, the high latitude caps predominate the world wide annual. 
variation, and from A/D<1, the low latitude belt predominates the world wide 
diurnal variation. “ 

From the above consideration it may be said that the upper conducting layer is 
not unique as hitherto considered, but somewhat complicated in structure. 

This paper was read at the meeting of the Society of Terrestrial Magnetism 
and Electricity of Japan held at Tokyo on May 1953. The writer wishes to express 

his thanks to Dr. Prof. Y. Tamura and Dr. H. Hatakeyama for their helpful encourage- 

ment. 
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Fluid Motions in a Sphere I 
Thermal Instability of a Rotating Fluid Sphere Heated within 


Ay 
By Tomikazu NAMIKAWA 


Institute of Politechnics, Osaka City University 
(Read Oct. 9, 1957; Received July 20, 1958) 


Abstract 


Thermal instability of an incopressible rotating fluid sphere with negligible 
viscosity heated within is examined. It is found that we are not able to have 


stationary convections ; Instability first sets in as overstability. 


J. Introduction 


In recent years the onset of thermal instability in horizontal layers of viscous 
fluid, heated from below under the influence of Coriolis acceleration and magnetic 
field has been investigated [1]-[4]. These studies have disclosed a number of novel 
features. Extent of these investigation to a fluid sphere was carried out by Chandra- 
sekhar [5] and to a rotating sphere by Takeuchi and Shimazu [6] and Chandrasekhar 
[7] in the case of axial symmetry under modified boundary conditions. We are — 
interested in the convection in the earth’s core or in the stars under the influence of 
Coriolis force and magnetic field. The present paper is devoted, as the first step, to 4g 


BO es of a rotating fluid sphere with mcersal heat sources. 


2. The Equation of the Problem 


a 
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where @ is the perturbations in temperature. 


We can elminate 6 from eqn. (3). by taking its curl. Making use of eqn. (4). 
we have 


a 


Ow du 
0x ui zi DOP ee: 
Gree tl + yo + aa (5) 
where w denotes the vorticity. Next, by taking the curl of eqn. (5), we obtain — 
~ 2 pu=reurl(poxr)—yp'at 205% (6) 


We can express the velocity field as a superposition of poloidal and toroidal vector in 
terms of two scalars. 
u=pT xi,+curl(pS xi,), (7) 
where i, is the unit vector of radious direction. Two scalars T and S are expressed 
as follows 
2= SL OY RandS= p> Sal tpt) awe, (8) 


where Y™ is a surface spherical harmonic. 
The following procedure is equivalent to take the poloidal and the toroidal 
components of eqn. (5). 
. Multiplying eqn. (5) by r, we get 


ore). =yp?(re,) +2.9r0# ae (9) 
and ‘ 
— SF ri(rur)= 776 —vy' (ru,) +22r—— oe, (10) 


? eft 4 


3 | We express at Zi8 On(r, t)Y™. ones vie 60: eqns. (2), (9) and (10) sade executing 


integration ri respect to 3 ons y, we ese 
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Ms and I2(7, =2T'Rlr, t), | (14) 


but for simplicity, we write 7’ and T7(7, ¢) as ry and T%(r, t). Thus eqns. (11), (12) 


and (13) becomes as follows. 


9 gm ct {|< 24 \_vins jon- one 
g ei mon d 23 
il 2 2 
n(n+1)a®?—_ —vn(n+1)— <r <—)- —n(n+1) 
— 2042{ = ream mM) yn-t ¢ S24) + n(n+2)(n+m+1) 
2 haa M 2n+3 
ieindae eee (16) 


and 


— n(n) vila ee <)- =n(n+1)| 22 2 Fontan 0] {= —( ae) 


—n(n1)} PS : —~20a? (a— a raf d TR 


£ n(n+2)\(n+m-+1) (n+2)_ 2 _ a bs . ye une) 
2n+3 a a Tae dr* r si 


=yn(n+1)a°o” ~ (17) 
Equations (15), (16) and(17) are the basic equations of this problem. 


3. Boundary Conditions = 


We treat eqns. (15), (16) and (17) neglecting the viscous term. This simplification” = 
esha ea we ose some geopiets or Sage ia PS For example, — 


scene at Heo ditto ow ei is ig 10-~ 10-4, which is + imc sm 5 
ete c of Coriolis. Yoree [9]. ye 
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Ti S3—T?} —...... (18) 


etc. 


We can satisfy the boundary conditions on 0%(r, ¢) and S%(r, t) by expanding 
them in terms of the various modes which form a complete set of orthogonal functions. 


Thus , 
fase?) CAD J 41(a57) 
o™ t Qn ahige ae ake d 2 Y, es ple n Ne 
ir, H= » ee aE an Se: uD Sa ar : (18) 
where /,,, 4 demote the Bessel functions of order Mii Zoek heres Ss (J=1,02, ca) are 


® itS zeros. 


5. An Examination of the Principle of the Exchagne of Stabilities 


Bisex 2 w3, -<-«| motion 

We shall first examine stability of [S%, T7%,,, ....] motion taking only the first 
two motions S? and 7%,,. It seems that this simplification does not alter the results 
to any great extent [6], [7]. We expand S%(7, t), T%.,(7, ft) and 6% (r, t) as follows, 


Sar, t) = Inst (a;7r) os Ine glean) 
fae. SSDS ey oe a, pins ee 
and 
: be I eae) 
ON(r, )= DOr. Der Saget | (20) 


i eee aj'S are the jth roots of Instlaj)= ii _ . 
ute te Tf all quantities very like e, eqns. (15), (16) and a7 are reduced to” 
Bscsena) feratls © Cat+aay\ORy—2Ban(n41)SB5=0, 2) 


DOME TE 200'{ —- — Niet La 2 54 +imT rs, ia (22), 


Bis ds ry Alastor gufizoienret “Withee Es Rat an Shgnne sit worl aw ATS 


om a . 
me " 


‘ania 
- 
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where 
Pi=ka;*, po= —2mkT, 
Es = n(n+2)(n+m+1)\(m+1—m) _ m jer? Aan Ranta. 2 (27) 
(1+1)?(2n+1)(2n+3) n(n+1)?(n+2) aj 
Substituting o=x+7y in eqn. (26), we get - 
+ pyx?—(3192+2Po.9—Qi + Ri)x—PiY’—Pi Poy+PiPi=0 (28) 
and mercy "6 
Le 
24H (2 asigt = pi, ait (qi) 1 4 ae ee 29 
(39+ po) x? + pil2y+ b2)x— YP —Po¥+(Qi—Ri)y Gali | @ ) 
As pi, @, a mT are positive, we have no steady state solutions. 


¢ (n+1)(m+2) 

i) axially symmetric solutions 

In this case, eqn. (26) is reduced to 

o+p,o? +(qi:-R)o+pigi=0 (30) 

When q,—;<0 this eqn. is able to have positive real roots, but when g,—R,>0 it is 
not possible to have positive real roots. Eqn. (30) allows a negative real root and the 
corresponding solutions lead to an exponential damping of an initial disturbance. 
We have to examine whether under these conditions the complex roots of (30) can 
have a positive real part. If x and y denote the real and the imaginary parts of 
complex roots of eqn. (30), they must satisfy the equations 


ee abt Pde Rie (31): ; 
and ae , 
x8 bes ih ee +P oie =0. Roo (32) a 
From the above two eqns. we have . ft ee 
| 84 8px? +20 Dita. Re—A.R, =0. Saupe 139) aa 


This eqn. admits a Positive real spot aaa Retna >), x we ‘shall have the ctital 
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When 


R 2(n+1)(n+2) | 
1> Ant in+2)—-1 ’ (38) 


we have postive value of x. 


When eqn. (26) has no positive roots, its comple roots must penehy 


Ob py — (92+ bey ~ai) +9(2y+ Po) + Ri) — Py? +poy—q1)=0 (39) 
and 
3y+ pox? + p,(2 — yy? ey eee ee ny 
Cr+ bast dQy+e:)e— Ho +Pey—g+Ry eI =0, (40), 
Egn. (39) admits positive roots when , 
yY+P2y—-Q>0 (41) 
Eqn. (40) admits a positive root when 
3y+p2>0, (42) 
and y 
ie ai R mreT = f 4 
WY + P2941) + bein »}<0 (43) 


General examination of the inequality of(41), (42) and S is troublesome, but | 
- these inequality can be clearly satisfied by the condition ms 
caesar dae a5 se (44) 
(2+1)(n+2) 
Thus we have a case of overstability when period of oscillation is smaller than 
a critical value. The critical number of vibrations is 
2m nt are ; pon eee (45) 
ow Dnt : : 
ete is smaller than the angular velocity of sphere f for ou value of x tends to 
zero when x tends to infinity. — Meee 
B) 1p Semen OR cine vis aT theme te amen oul 
follows nee ad 


y> 
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From above three eqn. we get 


ep a! gz (n+1)(n+2)? R| . 


(nln +1) o-imeT }{ (n+ I)(n+2)a3?—imas w+ ka; 


m(nt+2)(ntn+1(n+1—mM) , 2272=9 , (50) 


e (2n+1)(2n+3) 
Eqn. (50) is reduced to vee 
w+ (p,+ip2)w? +(q,— Ry ret AW ook Wie el: ce an Riaos (51) 


ov (nt+1)(m+2) p 
a," 
As eqn. (51) has the same form as eqn, (26), we can apply the examinations of 
eqn. (26) to eqn. (51). We have no steady state solutions. 


where ~;, f, and gq; are the same notations in (27); and R,’ 


i) axially symmetric solutions 

In this case the discussion of [S,, 7,4; ----] motion is entirely applicable; When 
R,>q,; we have an increasing amplitude instability, but when R,<q,, we get increasing 
amplitude oscillations and the critical pa ae of oscillations is equal to that of [S,,, 
Ta eee HOuUOns 

ii) asymmetric solutions 

The real roots of eqn. (51) must satisfy the equations 


: x84 pix? +(q,—Ry')x+h191=0, (52) 
and x 
24 pi pox + Tp =0, | 53 
Dox + Piet aig ea 1 (53) 


The discussions of [Sm, Tt, ....] motion is entirely applicable to [Tx Srai----] 
motions, letting n>n—1 _and R,>R,!.. When 


Ri'<q and ere ron (54) 
we are acase of overstability. The critical frequency of vibration is 
: Datos . ED OA ee 
er ss m1 we ae all 
and an cura, eit n tends to re . ae = bes 


oe ae ee ramerty aL <a, uses asl oe 
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7. Appendix | 


Equation of Motion [5], [7]. 

As we are seeking the effect of Coriolis force, we shall ignore the rotational 
flattening of the sphere in treating this problem. In the state of no convection, the 
temperature distribution, T(r), is governed by 


Ee ee (Al) 
ot 
where « is the coefficient of thermometric conductivity. The solution of eqn. (A.1) 


appropriate to the problem on hand is 
Ty= (a1) =8(0—-9). . (A2) 
Kk 


We have assumed that 7)=0 at r=a; This entails no loss of generality. 
If we assume the temperature distribution in the perturbed state 


T=T+6, (A3) 
the equation govering @ is 
0 (up) T= 46°0 ; ¥ (A4) 
where w denotes the velocity. 
» ou + o(up)u= =—pP+oyV+ovp*u+ ou x 2Q, (A5) 


where p, P, V and v denote density, pressure, gravitational potential and kinematic 
-viscisity respectively. 

Following Rayleigh [8], we shall allow in eqn. (A5) for ie variation of density 
only is so far as it modifies the external field. Thus in eqn. ts 5) we shall a 0 
whit occurs in front of 7V by a Si, a 

go =p aT), Soma AC) 
where a denotes the coefficient of volume expansion and Po is the density at. poe 
(where T=0) me pee 0 pegine elsewhere sap eles (A. 6) as po. wha m Te 

Now making use of rela beco: esomed Lata ay? 


. 


IMIS 


—_ 
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where 
ee 2h A wot 
and 
=-°9-y 4 Br (2atr? = ok (A12) 
00 


We shall suppose that w and @ are small quantities of the first order and that we can 
ignore products and squares of them. 
Then, the equations in the perturbed state are 
Ou _ 


a aoe +70r+yput2Ux 2 (A13) 


and 
20 = 670-+26(7U,), , (Al4) 
where Op is the perturbation in pressure. 
Boundary Conditions 


Boundary conditions which the solutions of eqns. (A.13) and (A.14) must satisfy 
depend on whether the bounding surface is rigid or free. but in all cases we requir that 


#=0- and 4v4,=0 (on a bounding spherical surface). (A.15) = 


Additional boundary conditions follow from the equation of SIDE and depend o on 
se the oe of the bounding surface. : 
_ Equation of continuity in spherical polar coordinates (7, 3, 9) is 


het 


Ou, u 1 Buy eal ae Tur go eainate d Saeed er 2 . 
Sete ee We OLY __1 _dOMy ayhh imeet ome ae 
ee Or rt 00g ae ying ao sank 
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Ge d Lie 
na+1)a* 8 —yn(n ti ya{ (PZ) nant] 2 
— 22a: A (a+ Dm) pn=1 : (St 1) ¢ 2AA(N+ M+) ,-cn42y 
2n— L PH) (A18) 
Z : x rs ms) | =0, 
| aS im and. 7=0 on rsl. 


dr 
(ii) on a free bouding surface y 
In this case we must have in addition to the condition (A.17), the viscous stresses 
P,»9 and P,» vanishes. 
The expression for these stresses are 


J PD ie daca 


Ou uU Ou 
P = v( : 7 aes °| 
PEt rsinvdg  r Or / 
Since w, vanishes on a bounding surface, the vanishing of P,» andP,, requires 


Dug _ Uy _ Op _ Ue (A.20) 
Or r Or Yr 


According to eqn. (7) these conditions are equivalent to 


(A19) 


7 (2 or) witds, & 1/83 — 2 oo) = 0 

GU Fe ar rar ¢\ ar Yr ar 

and ; 
dln 0 on r=l.- > (A21) 
dr 


Since S*¥=0 on 7=1, the first. of (A. 21) can be expressed BS ee coat 
ae | mate Bay eee oe hee 
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Thermal Instability of a Conducting Fluid Sphere Heated 
within under a Uniform Magnetic Field 


By Tomikazu NAMIKAWA 


Institute of Politechnics, Osaka City University 
(Read Oct. 9, 1957; Received July 20, 1958) 


Abstract 


In this paper the problem of thermal instability of non-viscous fluid sphere 
heated within under a uniform magnetic field is discussed. Differing from the 
former case of thermal instability of a rotating sphere, there exist axially 
symmetric steady state solutions and only (U,”", V’ns;, ..».) motions can arise. 
Calculation suggests a relation between perturbations in magnetic field and vorticity. 
Under normal terrestrial conditions, instability arises as ordinary cellular 
convections, but under astrophysical conditions, it can arise depending on intensity 
of the uniform magnetic field, as cellular convections or as overstability, oscilla- 


tions of increasing amplitude. 


1. Introduction 


In paper [I], we have examined thermal instability of a rotating fluid sphere 
and found that there is no steady state solutions; instability first sets in as over- 
stability. We shall examine in this paper the influence of magnetic field on thermal 
convection in a sphere. 


2. Equation of Motion 


We shall consider a homogeneous conducting fluid sphere of radious @ under a 
uniform magnetic field, the direction of which is taken to be z- axis. Assumption 
about the source of heat generation is the same as that of paper [I]. 

The equation of motion, of Maxwell, and of heat conduction appropriate to the 
problem on hand are 


p H+ o(u-pyu=—PP+ ous or V+ nd XH, (1) 
curl H=47J, (2) 
curl#= — OH (3) 

ot 
div H=0, div E=4nc’q, (4) 
J=0(F+yu~x HA) (>) 


(193 ) 
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and 


A+ (wp)T=op'T, (6) 


where E, J, H, o and » denote the intensity of electric field, current density, intensity 
of magnetic field, electrical conductivity and the magnetic permeability. The meaning 
of other notations is equal to those in paper (ll. 

The equations governing small departures from the stationary state are (Appendix 


of [1]; [2]) 


99 499 4 26(rur), (7) 
ot 

ou il uw oh 

eee Sey Pp 2 

Ar = Or “ae (8) 
pha Glee 
= 1 soph (9) 

and 

divu=divh=0, | (10) 


where H and A denote magnitude of uniform magnetic field and perturbation in 
magnetic field and 


= . eet 7 ‘ _— 
| Vey = (11) 


Taking the curl of eqns. (8) and (9), we have 


ars ee : kg 
3 ye BAP Be gir oe dz TER Me (12) 
and — 
= ae s eae tea ~ Vine sr _Jji)- seqeue) = 
: 3 “yy 0 = sre) tad: bean ar pee 


Ss ifs Aan, aking the cof ea. 12, we ast Sere ids haan ate oN asa 


Fluid Motions in a Sphere II 195 


Multiplying eqns. (9), (12), (13) and (14) by 7, we get 


Orhr) — Fp OW + apr rny), (18) 
ve O02 
9 (rer) = alt, £ ay (curk h)+vp(ro,), | (19) 
ot “Arp 
eOierr hi) = Soh ames h),., (20) 
Ot Oz 
and 
~ 2p (ru) =7L9— HF pop _yye(ru,), (21) 
ot 409 Oz : 


Operating Y% to eqns, (7), (18), (19), (20) and (21) and expecuting the integration 
with respect to @ and ¢, we get 


Ss Le yf (n=) —n(n+ js 


(n— Fens m) nm-1 d —n[ fm n(nt+2)\(ntm+1) .-¢ (+2) 
: ~ nf =D o df ae ne eas 


: te8 
x SU) = ue 0 ? (23) 


ee —yn(n+ p2{2( nd.) —n(nt 1} ue 


_ Ha (n— mire m) fre ad agit Tm AC ant2)(n +m-+ 1) p(s) 
rod at 2n-1 2n+3 
sae aly Gn a ? 3 (24) 
a es Be TS ind ee ee se 


= f(A -nne 2 
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3. Bouudary Conditions 


As in [1], we shall discuss eqns. (24) and (25) neglecting the viscous term. 


We must require that 
ee =UsD=0, (27) 


Additional boundary conditions are the continuity of magnetic field and tangential 
components of electric field on the boundary surface. 


4. Solutions of Basic Equations 


From eqns, (22), (23), (24), (25) and (26), we find BREN there exist groups of fluid 
motions and magnetic field as follows 
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We shall examine the (U™, T™, V%.,, S™.,; ....) motions taking only the first 
four terms in considerations. Simple calculation shows that there is no (V™, U™,,, 
S”, 7T™,,, ....) motion because continuity of the magnetic field cannot be satisfied on 
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E=y curl H—pu~ H. (33) 
Making use of eqn. (29), the internal electric field is calculated from eqn, (33): 
The external electric field must satisfy eqns. 
divE’=0 and — 1 OF curl E, (34) 
and exprssed as follows: . 
SH Bar YY 2) 4 Czp(r- wy), . (35) 


From egns. (33), (34) and (35), it is found that the continuity of tangential components 
of electric field is satisfied. 


If all quantities vary as e*’, substituting relations (29), eqns. (22), (23), (24), (25) 
and (26) are reduced to 


(Aa? +0,;2)0™ ;—2Ban(n+1)U™;=0,  - (36) 
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symmetric solutions of marginal stability, [U,;°, V2", 7°, S2°, ----] motion, letting m>0 
in (43). 
Thus Roos (n+2)a;'Q : 
6 n(2n+1)(2n+3)(n-+1) 
— 95° Loney rn 
iMag : (44) 


6. An Examination of the Principle of the Exchange of Stabilities 


Here, we examine axially symmetric case of eqn. (41). 
In this case eqn. (41) is reduced to 
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Eqn. (45) has no real positive root, but eqn. (46) has real positive root if 
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If : 
bq-r<o, (55) 
eqn. (54) admits a positive real root. 
Thus if . 
q>0, rs0 and pq—r<0, (56) 


the principle of the exchange of stabilities can not be applied and we have a case of 
overstability. 


According to eqn. (51), the condition (56) ds equivalent to 
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’ Rearranging the outer inequalities of in or (58), we get 
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7. Onset of Instahility in a Magnetic Field under 
Astrophysical Conditions 


a 


The inequality 
e>n oe: (66) 
is characteristic of astrophysical conditions [2]. When (66) holds, the inequality (61) 
or (62) can be simplified to 


Q> Ont VEnt) 4 2=Q, : (67) 


From eqn, (65), we can write for a considerable rarnge of Q satisfying (67), 
6 2 
Roa ae ( ees ). 68 
ae OS Ee (68) 
This is independent of the strength of the prevailing magnetic field. In axially 
symmetric case m=1 and minimum value of a; is a,: thus when (66) holds, ae 
in the form of oscillation of increasing amplitude can arise for 


) Q>54a," . : (69) 
and the critical value R®°, remains practically constant and has the value 
— GP Ke 
rae | Q< 55%) ; +0 (70) 
For Q> “54,2, R°®, can be written 
2 at Fate, >: 
7 \2 
ied, = con (2) P. 
t Kee \ =a (71) 
where ; é 4 
Reo Pras @ w ‘ 
at y 10 se : bg ; ; aj : 
= r 3 os et y, yf i ae - 7 "77 2 7 >} 
On the other hand, if * han Lhe See 2 a oe ‘jagiat, eg aes 


sil Q<5a;%, "2 5 Lee 
ge ee te ee “inesgatt alt ike Soaiteupane 


the principle of the exchange of stabil lite ss will be a le and instability s 


Fluid Motions in’a Sphere IT 201 


From eqns. (11), (42) and (75) 


2 
Q=V2e. (76) 
Me 0 
Thus we find 
peels 77 
it (77) 
where oF. is the time required for a magneto-hydrodynamic wave to travel a distance 
equal to 2rvJ/5 g. , 
ay 
Letting 
2 
Q=sai= Vi (78) 
we get 
=> V5 a7 . ; ; (79) 


_ Thus we may summerize that so leng as the velocity of magneto-hydrodynamic wave 
is less than V,, instability will arise through cellular convection, but for 


Vi<V< Vou—, we will have overstability when R° reaches the value giyen by (70). 


8. Examinations of Onset of Instability in the General Case 
Using eqns. (65) and (72) we can draw in the (ie Q)— plane the cree value 


of R. The convection curve starts on the R- axis at the pons a an 
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value of Q at the point where the two curves intersect is given by 


(fie: ST ae Q, 
Q, ae es (86) 
Kk Kk 


Thus for Q<Q,’, instability arises as a cellular convection, but for Q>Q',, we shall have 
overstability when the critical value of R° is reached. 


9. Concluding Remarks 


It has been found that only (U%, V4, 7%, Stu, .-.-) motion can arise in a 
non- viscous sphere under a uniform magnetic field. The remarkable feature is that 
the order and degree of the spherical surface harmonics of perturbations in magnetic field 
and vorticity are equal. This may suggest a relation or resemblance between the 
magnetic field and vorticity, which is first proposed by Batchelor [3]. 

Under the normal terrestrial conditions (y7>«), the inhibition of convection by a 
magnetic field is a very pronounced effect, but under astrophysical conditions («>~7), 
this is not. The results are quantitatively ‘equal to those of Chandrasekhar’s investi- 
gation of instability of a layer of fluid heated below under a uniform magnetic field. 
We shall treat the problem of instability of a fluid sphere under simultaneous action 
of rotation and magnetic field in the next paper. 
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Thermal Instability of’a Rotating Fluid Sphere Heated Within 
under a Uniform Magnetic Field 


By Tomikazu NAMIKAWA 


Institute of Politecnics, Osaka City University 
(Read Oct. 9, 1957; Received July 20, 1958) 


Abstract 


The influence of simultaneous action of Corios force and magnetic field on 
convection in a non-viscus fluid sphere heated within is examined. It is found 
that there is only axially symmetric solutions of the marginal stability, and 
inhibition of convection by a magnetic field is pronounced when the intensity of 


magnetic field is larger than a critical value. 


1. Introduction 


In paper [1] and [2], we examined the influence of rotation and magnetic field on 
the thermal instability of fluid sphere heated within. We are interested in the origin 
of the earth’s magnetic field and stellar magnetic field. According to the Elaasser- 
Bullard’s dynamo theory of cosmic magnetic field, the maintenance of a magnetic field 
is induction effect in the earth’s core or in the stars in which a thermal convection 
may exist. As the earth and the stars are rotating, we must solve the equation of 
fluid motion, of Maxwell and of heat conduction, taking into consideration Coriolis 
force. 

As this is a very difficult problem, we treat in this paper, these equations in 
linearized form and examine the influence of a simultaneous action of magnetic field 
and rotation on the thermal instability of a fluid sphere heated within. 


2. Equation of Motion 


Consider a homogeneous conducting fluid sphere of radious @ rotating with an 
angular velocity about z-axis under a uniform magnetic field, the direction of which 
is taken to be z-axis. Assumption about the source of heat generation is the same as 
that of paper [1] and [2]. The equation of motion, of Maxwell, and of heat conduction 


appropriate to the problem on hand are 


# - +o(u-p)u=—Pp+20ux 2+ oyp'u+ op V+ ud x A, ee 


curl H=42J, (2) 
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cit H=—y ae (3) 
div H=0, (4) 
div E=4nc’q, (5) 
J=o(F+pnux HA), . (6) 
and of +(u- pT =«p°T. (7) 


The meaning of notations is equal to those in paper [1] and [2]. 
The equations governing small departures from the stationary state are ((1], [2], [3]) 
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Taking the curl of eqns. (9) and (10), we get 
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Operating Y” to eqns. (8), (15), (16), (17) and (18) and executing the integration 
with respect to @ and ¢, we have 
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The boundary conditions are the same as that of [2]. The vanishing of perturbation 
in temperature and radious component of velocity, and continuity of magnetic field 
and tangential components of electric field on the boundary surface are required. 


4. Solutions of Basic Equations 


As in [2], we shall examine the (U%, rT, v™, S™,,..) motions taking only the first 
four terms in consideration. There is no (V™,S™”,U%,,, T%4,..) motion because the 
continuity of the magnetic field cannot be satisfied on boundary surface. 

Expansions of velocity and magnetic field inside sphere or outside sphere, are 
equal to eqns. (29)-(35) of [2]. 

If all quatities mary like e*’, basic equations are reduced to 
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and 
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and 
R/Q=| aekae ee ee, (a? as (40) 
n+1 (2n+1)(2n+3) n* 2(m +1) 
As the right hand side of eqn. (40) is negative for all value of (=m), we have 
no steady state solutions in general case. 


Letting m—>0 in (39), we get 


ee, ee fs Nyay n(n +2)2(n+1)? I= (41) 
n(n+1)frb (n+1)? (2n+1)(2n+3) Q 

We have only axially 
symmetric solutions of marginal | t=10° 
stability, the (Ut, V2- 7°,S3,..) i2b 
motion. T=10° 

The minimum value of alt 
aj; 1S a,=4.5; critical value of R T=103 
is given by [ 

Rom =a,4 ere (42) 


Re has its minimum value 


when 


Q=5./2T. (43) 


6. Concluding Remarks 


0 | 2 3 4 5 fog Q 
It is found, from eqn. (42) Dependence of the characteristic number FR for the 
that when angular velocity and onset of convection on the number J and Q. 


other constant of fluid, x, 7, 0, etc., are given, the inhibitation of convection by a 
magnetic field is pronounced only when intensity of magnetic field is larger than a 
critical value. The behaviour of R resembles to that of Chandrasekhar’s investigation 
of instability of a layer of conducting fluid under simultaneous action of Coriolis force 
and uniform magnetic field. 

We shall report in the near-future the results of the examination of non-steady 
state. It is supposed from the results of this paper that it is very interesting to 
investigate the convection in the earth’s core, solving equations of motion, of Maxwell 
and of heat conduction in more complete way. 
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LETTERS TO THE EDITORS 


The Result of Observation of the Rate of Ion Pair Production in 
the Atmosphere during the Solar Eclipse, Apr. 19, 1958. 


The observation of the rate of ion pair production in the atmosphere near the 
ground has been carried out by the auther at Tanashi, Tokyo, from the view point of 
the atmospheric physics. 

The instruments used for this observation are the ionization chamber and the 
vibrating reed electrometer. The form of the ionization chamber is a cylinder 90cm 
long and 20cm in radius having a volume of 27 liters. The wall of the chamber is 
myler of 7 micron thick (8.5x10-4g/cm?). The inner electrode, 5 mm in diameter, 
is supported centrally in the base of the chamber by a shielded insulator (teflon), and 
is connected to a vibrating reed electrometer (made by “ Applied Physics” U.S.A.) by 
a shielded wire. The vibrating reed electrometer measures the saturated ionization 
current, to indicate the rate of ion pair production. The ionization chamber is set in 


the air flow cylinder. The air intake is placed at 1m height from the ground surface. - 


19, 4, 1958. (Tanashi, ToKYo) 
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The detailed discussion on the result of observation mentioned above will be 
reported in the future. 


By. M. KAWANO 

3 Electrotechnical Laboratory, Tokyo, Japan. 
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On the Magnetic Field of S, in the Middle and Lower 
Latitudes during the II Polar Year 


1. Summary 


The geomagnetic S, field was calculated as the mean of five disturbed days in 
every month during the II polar year, by the same method which was adopted on the 
analysis of Sy [1], using the data of 45 observatories in the middle and lower latitudes. 

Comparing this result with that of S, calculated by H. Maeda, [2] Sp field was 
computed. 


2. Brief Description of the Method 


Using all available data, values of X components at the cross-points of parallel 
circles every 15° in latitudes and meridians every 30° in longitudes, the treatment was 
carried out about every two hours, i, e, 0%, 2”,--+++ , 22" LT, in order to obtain 
values in the mean state. Potential of Sy, field of 5 (S+ W) was calculated through a 
method of graphical integration using the above X components. Values in latitude 
60°N were nearly equal to values in latitude 60°N obtained from the analysis in polar — 
region [3]. 


Potential V is expressed in the next equation. 
V=ry3P72(0) {A% cos (T+4)+B* sin (T+A)}. 


ry: radius of the earth, T: universal time, 4: geomagnetic longitude T+A=t¢: local 
time, P7.(0): normalized function of Ad. Schmidt. 
The value of A”,B® are indicated in Table I, In table II are indicated those of 
_ Sy obtained by H. Maeda [2]. aS 
Comparing these two Tables, it is Cate 
equal and the comparatively large te: ae 
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Taking the difference between the coefficients of S; and S,, we find the values 
of Aj, A?, Ai, A3, are very small and those of Aj, A, ;B3, Bi, Bi, Bi, B} are comparatively 
large. 

According to this table, it is seen that the diurnal and the semidiurnal terms 
predominate. From the values in Table III, the map of equipotential lines are drawn, 
which coincide nearly with that of S,*obtained by S. Chapman. 


Table II. The Coefficients of Sg. 
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Fig. I. Equipotentials of the Sg field, where @ shows position of the Sun. 
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Fig. II. Equipotentials of the Sg field, where ne shows position of the Sun. 
[after H. Maeda] 
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Prevailing Wind in the Ionosphere and Geomagnetic Sq Variations 


Recent observation of the wind motion in the ionosphere [1] has detected the 
considerably great velocity of the prevailing wind in the E region of the ionosphere. 
The prevailing wind in the ionosphere, as the zonal wind in the troposphere, is 
predominantly in the east-west direction and constant in local time. 

The author derived the horizontal wind system in the E region from the 
geomagnetic S, variations [2] [3]. In his calculation it was assumed that the wind 
motion corresponding to the S, variation is periodic with local time. The prevailing 
wind motion, constant in local time, might make some contribution to the S, 
current system due to the daily variations of the ionospheric conductivity. If the 
effect is noticeable, the treatment of the author would be based on the irrelevant 
assumptions and the results should be revised. The purpose of the present report is to 
show that the prevailing wind would play no role in the first-approximate treatment by 
which the previous calculations were carried out. In the first-approximation we assume 
that 1) the S, current flows horizontally in the £E region, which is a thin spherical 
shell, and 2) the geomagnetic field is that of a centered dipole alined along the 
geographical axis. Further the prevailing wind is assumed to be axisymmetric as the 
zonal wind in the troposphere [4]. Namely, the wind velocity v is 


u=0 , 
imme 
V=WG) 5 


where u is the southwards component, v the eastwards component of v and @ the 
colatitude. The eastwards component v is, in general, a function of 6. The electric field 
E is composed of the induced field and the static field as 


E=voxH+grad S, (2) 
where H is the geomagnetic field, S the static potential. Ohm’s law is given as 
E=(r|J, (3) 


where J is the current density and [7] the resistivity tensor defined as the reciprocal 


of the conductivity tensor [o] viz. 


[a] = Oxx Oxy [7] = [o}"* a Vx Vay 
(4) 
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AE ae Op Mar cee ear ET 
2 
Oxx yy t Oxy 


(215 ) 
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In (4) an orthogonal x, y, and z coordinates are taken to be southwards, eastwards, and 
vertically upwards respectively. Next we transform the coordinates so that the « axis 
coincides withe the direction of J. Then, by denoting the transformed quantities by a prime 


As cosa ae [7] aa i 


—sing cosa sina cosa (5) 

hes be costa t Fyy sin’a | fi (6) 
, 

E'y= [ —(n—Tyy) Sina Ccosa— rey) i oe (7) 


where « is the angle between the x and x’ axis. This angle is, of course, a function of 
6 and the longitude in general. 


Now we take a line integration about a closed current path, provided that 


divJ=0. : (8) 
Then, 
f Has= § [r\Jds = -tgtDh 
current current 4 : j - 
where § is the path of the electric current. By (2) and Stokes’ theorem abate A 
¢ Eds= f (wx H)ds= ff rot (OXE)dS.. os » Da 
current current f : = % 


where the double integral means the integral over any regular surface boutidied by s. - 


By (1) | 4 7 
te ) ry cL : cette at = = aie E 


7 , : - ré IE. > Je = — = (6 f 
teh Sos 2h is ce ec erent Oak sae | 


——— 2s 
= Kine = i 


(Bosse 
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Jee (15) 
while 
Len COS’a+Tyy Sin*a>0 . (16) 
Therefore 
Z =O (17) 


From the above argument it is concluded that the axisymmetric wind motion 
produces no electric current under any global distribution of the ionospheric conductivity 
in the dynamo theory, based on the above two assumptions, which have usually been 
made. The induced field is cancelled out by the resulting static field. Otherwise, 
free charges would accumulate constantly around the poles. 

Actually, the two assumptions 1) and 2) would not exactly be met. The depen- 
dency of H. on the longitude leads to production of a current by the axisymmetric wind 
motion. The prevailing wind actually contains a small southwards component which may 
give rise to a current flow. This has not been considered so far. Further, the 3-dimensio- 
nal treatment will yield a different situation. However, these possibilities are likely to 
arise in the more detalied treatment, say, in the second order approximation. The 
outstanding character of the wind system corresponding to the S, variations is 
represented in the results of our calculations [2] [3] [5] [6]. 
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Meeting of the Society of Terrestrial Magnetism and Electricity : 


The 23rd General Meeting was held at the Tokyo College of Science on may 17-19, 
1958. 


Number of the Reports read at the Meeting: 
Rock Magnetism, 12; Atmospheric Electricity and Atmospherics, 8; 
Cosmic Rays, 10; Ionosphere, 11; Geomagnetism, 16; Earth Current, 2; 
Night Airglow, 3; Special Reports on ‘Antarctic Research Expedition, 3. 
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